function Shell_mixed_FEdemoshrt_grID_wtp_tomlab_varagin_j_bs(subsample,inc,set_seed,matlab_seed,get_nb_draws)

%matlab_seed=1353;
rng(matlab_seed)

	tmp_file1=strcat('\\c3\rdat\SHoude\Research/sears/estar_data\refrigerators\choiceset_identifier_trimester_week_zipcode_2008_2012_v11022017_struct_',num2str(subsample),'_',num2str(inc),'_seed_',num2str(set_seed),'.csv');
    id=load(tmp_file1); 
	Nmax=size(id,1);
	bs_max=10*ceil(Nmax/20000);

parfor j=1:bs_max	
	%index_bs=[(1+(j-1)*(1000)):min(j*(1000),Nmax)];	
	index_bs=[(1+(j-1)*(2000)):min(j*(2000),Nmax)];
	Master_hom_FEdemoshrt_grID_wtp_tomlab_varagin_j_bs(subsample, inc, set_seed, j, index_bs);
	Master_mixed_FEdemoshrt_grID_wtp_tomlab_varagin_j_bs(subsample,inc,set_seed,matlab_seed,get_nb_draws, j, index_bs);
end

end